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ABSTRACT 

The  fracture  problem  of  laminated  plates  which  consist 
of  bonded  orthotropic  layers  is  studied.  It  is  assumed 
that  the  medium  contains  periodic  cracks  normal  to  the  bi¬ 
material  interfaces  and  the  external  loads  are  applied  away 
from  the  crack  region.  The  field  equations  for  an  elastic 
orthotropic  body  are  transformed  to  give  the  displacement 
and  stress  expressions  for  each  layer  or  strip.  The  un¬ 
known  functions  in  these  expressions  are  found  by  satisfy¬ 
ing  the  remaining  boundary  and  continuity  conditions.  A 
system  of  singular  integral  equations  is  obtained  from  the 
mixed  boundary  conditions.  Three  cases  are  considered: 

a)  The  case  Of  internal  cracks 

b)  The  case  of  broken  laminates 

c)  The  case  of  a  crack  crossing  the  interface. 

The  singular  behavior  around  the  crack  tip  and  at  the  bi¬ 
material  interface  is  studied.  It  is  s’lown  that  the  crack 
surface  displacement  derivative  has  a  power  singularity  for 
practical  orthotropic  materials  when  the  crack  touches  the 
interface,  i.e.,  for  case  (b).  In  studying  the  singular 
behavior  at  the  bimaterial  interfaces  in  case  (c),  it  is 
found  that  for  some  orthotropic  material  combinations  there 
is  no  singularity  in  the  crack  surface  displacement  deriva¬ 
tives  and  the  stresses.  In  each  case  the  stress  intensity 
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factors  are  computed  for  various  material  combinations 
and  various  crack  geometries.  The  results  for  orthotropic 
materials  are  discussed  and  are  compared  with  those  for 
Isotropic  materials. 
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1 .  INTRODUCTION 


In  structural  design  one  of  the  most  important  consid 
erations  is  the  fracture  of  individual  components.  Al¬ 
though,  fracture  may  not  always  mean  total  failure,  it  is 
considered  in  modern  engineering  as  an  important  problem 
for  safe  and  economic  design  of  structures.  It  would  be 
very  attractive  to  develop  special  types  of  designs  for 
which  the  structural  resistance  to  fatigue  crack  propaga¬ 
tion  is  improved.  In  the  aerospace  industry,  the  use  of 
composite  sheet  materials  with  buffer  strips  parallel  to 
the  main  load-carrying  laminates  seems  to  be  such  a  design 
practice.  The  process  of  manufacturing  composites  gives 
the  opportunity  to  improve  the  structural  resistance  to 
fatigue  crack  propagation  by  strengthening  the  material  in 
certain  directions.  The  increasing  use  of  composites  in 
structures  generates  new  problems  for  the  structural  de-  . 
signer.  Among  these  problems,  we  are  mainly  interested  in 
the  fracture  of  layered  composite  materials. 

There  are  two  main  problems  in  studying  the  fracture 
of  composites:  the  development  of  an  appropriate  failure 
criterion  and  a  mathematical  model  for  the  calculation  of 
the  related  load  factor.  The  failure  criterion  affects 
the  course  of  the  analytical  work  in  the  sense  that  it  is 
the  failure  criterion  which  generally  determines  the  phys¬ 
ical  quantities  that  one  should  compute  (such  as  the  stress 
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Intensity  factor,  the  strain  energy  release  rate,  COO, 
etc.).  There  are  many  failure  criteria  or  theories  which 
are  used  to  predict  failure  of  structures.  In  Elastic 
Fracture  Mechanics  where  only  small  scale  yielding  is  al¬ 
lowed,  K  <  Kjq  is  such  a  criterion.  In  this  case  failure 
occurs  when  the  calculated  value  of  the  stress  intensity 
factor  reaches  a  critical  value,  Kjq,  which  can  be  deter¬ 
mined  experimentally  as  a  material  property.  There  are 
also  other  one-parameter  failure  criteria  (such  as  critical 

plastic  stress  intensity  factor  K  and  J  integral)  which 

pc 

have  been  recently  proposed  to  predict  failure  from  elas¬ 
tic  to  fully  plastic  range.  K  is  a  very  highly  effective 
correlation  parameter  in  studying  the  fatigue  crack  prop¬ 
agation  phenomena.  In  aerospace  structures  the  basic 
problem  is  the  nucleation  and  propagation  of  fatigue  crack 
which  may  eventually  reach  a  critical  size  causing  cata¬ 
strophic  failure.  That  is  why,  in  this  study  we  focus  our 
interest  to  the  computation  of  the  stress  intensity  factors 
and  in  the  investigation  of  the  singular  behavior  of  the 
stress  state  around  the  crack  tips. 

In  studying  the  fracture  problem  of  composites,  a 
mathematical  model,  which  will  reflect  the  geometrical  and 
physical  properties  of  the  medium  and  the  real  mechanism 
of  fracture,  is  needed.  Because  of  mathematical  diffi¬ 
culties  and  the  lengthy  computation  that  the  analysis 
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requires.  In  the  recent  studies  the  geometry  and  the  ma¬ 
terial  properties  have  been  considerably  simplified.  The 
problem  of  a  multi-layered  isotropic  medium,  which  consists 
of  many  layers  and  where  a  crack  normal  to  the  interface 
can  appear,  has  been  treated  by  Hilton  and  Sih  [1].  In 
this  problem  the  geometry  is  simplified  to  a  single  layer 
between  two  dissimilar  half-planes  where  the  elastic  prop¬ 
erties  are  averaged.  The  same  problem  has  been  considered 
by  Bogy  [2].  The  problem  of  a  broken  laminate  between  two 
half  planes  has  been  investigated  by  Ashbaugh  [3]  and  Gupta 
[4].  The  extension  of  the  problem  treated  by  Hilton  and 
Sih  to  orthotropic  media  has  been  solved  by  Arin  [5].  The 
fracture  problem  of  a  composite  plate  which  consists  of 
parallel  load-carrying  laminates  and  buffer  strips  has  re¬ 
cently  been  solved  by  Erdogan  and  Bakloglu  [6].  In  this 
work  the  load  carrying  laminates  and  buffers  are  consid¬ 
ered  to  be  isotropic  and  linearly  elastic.  The  orthotropic 
case  of  the  problem  treated  In  [3]  and  [4]  has  also  been 
solved  by  Arin  [7]. 

The  objective  of  this  work  Is  to  Investigate  the 
fracture  problem  of  composite  plates  containing  periodic 
buffer  strips.  The  laminates  and  buffer  strips  are  assumed 
to  be  linearly  elastic  and  orthotropic.  In  general,  this 
Is  the  case  in  the  actual  plate  and  shell  structures  such 
as  those,  for  example,  which  consist  of  boron-epoxy 


composites.  It  is  also  assumed  that  the  fatigue  cracks 
may  appear  and  propagate  in  main  laminates,  in  buffer 
strips  or  in  both  normal  to  the  interfaces.  The  external 
load  is  applied  to  the  plate  parallel  to  the  strips  and 
away  from  the  crack  region.  Three  different  problems  are 
studied:  the  internal  crack  problem,  the  case  of  broken 
laminates  and  the  case  of  a  crack  crossing  the  interface. 
A  general  formulation  of  the  problem  is  given  for  plane 
strain  and  generalized  plane  stress  cases  by  the  use  of 
Fourier  Integral  Transform  Technique.  The  singular  be¬ 
havior  around  ends  and  at  the  bimaterial  interfaces  is  ' 
studied.  The  resulting  singular  integral  equations  are 
solved  numerically  and  the  stress  intensity  factors  are 
calculated  for  various  crack  geometries  and  various  ma¬ 
terial  combinations. 
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2.  ELASTICITY  OF.  M  ANISOTROPIC  ELASTIC  BODY 

For  an  anisotropic  elastic  body,  in  the  absence  of 
body  forces,  the  equations  which  relate  the  field  quanti¬ 
ties  can  be  written  as  follows: 

2.1  The  Equi 1 i br i um  Equations 


!l21  + 

+ 

3x 

3y 

3Z 

3x 

ay 

3Z 

3^z  ^ 

atyz 

+ 

3x  * 

ay 

3Z 

2.2  Strain-Displacement  Rel ations 

6,1“  e=5v  e=>^ 

3x  ’  '^y  3y  ’  3z 

Y  :.lV  +  lW  Y  ,111  +  1“  Y  =l“  +  ll 

'yz  3z  3y  ’  'xz  3x  3z  ’  'xy  3y  3x 

(2.2) 


2.3  Stress-Strai n  Relations  J 


p  ^ 

P  - 

•^X 

^y 

i.j  =  1.6 

°z 

-  fa  1 

^z 

^yz 

'bu 

^z 

* 

A . .  =  A . .  r 

^z 

^xz 

U  J1  E 

^xy 

^xy 

(2.3) 
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2.4  The  Field 


for  an  Orthotropic  Body 


For  an  orthotropic  solid  the  matrix  [Ajj]  is: 


[*u]  • 


^11 

Ai2 

*13 

0  0 

^12 

CM 

CM 

*23 

0  0 

‘l3 

*23 

*33 

0  0 

0 

0 

0 

*44  ° 

0 

0 

0 

°  *55 

0 

0 

0 

0  0 

Defining  the  inverse  of  [A.-]  by 

■  J 

[a..]  =■  [Aij]-* 


for  orthotropic  materials  we  have: 


(2.4) 


E  V  =  E  V 
X  yx  y  xy 


E  V  =  E  V 
'•y  zy  z  yz 


^z'^xz  “  ^x^’zx 


(2.5) 
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Substituting  (2.2)  into  (2.3)  and  using  (2.1)  and  (.2.4), 
the  stresses  and  the  equilibrium  equations  can  be  expressed 
in  terms  of  the  displacements  as  follows: 


0  34  ^+A  ^ 

®x  ''ll  3x  *  ''12  ay  ''13  az 


o  3A  ^+A  i^+A  ^ 
®y  ''12  ax  ''22  ay  ''23  az 

a  =■  A  M  +  a  ^+A  ^ 
®z  ''13  ax  ''23  ay  ''33  az 

^yz  *44' az  ay^ 

^xz  *55'ax  az^ 

^xy  *66' ay  ax^ 


(2.6) 


3*u 


+  A  ^+A  l!“+iA  +Al  +  fA  +A  1  =  0 


"11  ^  ^  "66  aF  "55  ^  ^  '"12^"66'  W  ^  '"13^"55'  3)^ 

a*«  3*V  ,  a*  _ 

ar^  ''44  aP  ■  '"i2'”66'  axay 


A  3!v  +  A.,liv  +  A  1!X+(A+A1  +  fA  +A  )  =  0 

66  ^  ''22  aF  44  aF  '''12  ''66'  3xay  '''23  ''44'  ayaz  ” 


a*w 


+  A 


3*w 


+  A 


a^w 


+  (A,,+Anc) 


a^v  _ 


'55  ^  "44  ^  "33  ^  '"13  ”55'  axaz  '"23  "44'  ayaz 


{2.7} 


2.4.1  Case  of  Plane  Strain 


For  the  plane  strain  case  we  have: 

u  =  u(x,y)  ,  V  =  v(x,y)  ,  w  =  0 
and  from  (2.2), 


(2.8) 
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e-.lH.  c-lX  V  -M+iZ 

3x  *  ~  ay  ’  'xy  '  ay  3x 


E  *  0  •  Y  “0 

"  •  ’yz 


Y  =0 
'xz 


(2.9) 


Thus,  thz  stress-displacement  relations  and  the  equilib¬ 
rium  equations  become: 


^  .  3u  .  ,  3v 

®x  *11  71  ^  *12  ay 

®y  *12  3x  *22  ay 

’“  +  A  ^ 


®z  * 

*13 

au 

3x 

It 

>» 

K 

*66 

/  3u 
'ay 

3 

T 

=  0 

yz 

XZ 

(2.10) 


A  ^  +  A  liH  +  fA  +  A  1  gjy  »  0 
''ll  IJF  ''66  ap  '''l2  ''66'  axay  ^ 


a*v 


+  A. 


a^v 


+  (A,-  +  Ace) 


a^u 


'66  ax^  "  "22  ^  '  '"12  "66'  3xay  ’  ' 

2.4.2  Case  of  Generalized  Plane  Stress 


(2.11) 


In  this  case  since  =  0  from  (2.3),  for 

the  average  stresses  and  strains  we  can  write: 


'xy 


[5] 


xy 
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where 


Ci]= 


“  r^  r4: 


X  -  -X 


and  A  - 


-  ''xyV^- 


(2.12) 


The  equilibrium  equations  reduce  to: 


^11  llF^A33|y^MA,2-^A33)^=0 

A  .5 _ X-  +  A  _ X  +  iA  +  A  ^  ^ ^ U  _  ^ 

''33  JyJ  ''22  '''l2  *33'  3x3y  '  ° 


(2.13) 


Considering  the  structure  of  equations  (2.11)  and  (2.13), 
the  equilibrium  equations  can  be  written  for  plane  strain 
and  generalized  plane  stress  cases  in  the  following  form: 


a  i-Ji.  +  i!“  +  a  -  n 

“1  Jx^  3y3x  ■  “ 

ifV  ^  g  ifV  g  3^U  _  g 

33F  3x3y  '' 


(2.14) 


where 


g  ,  "11 


6  - 


B3  =  1  +  —  for  plane  strain 
66 


B  = 

^  A33 


a  =  ^ 
^  A33 


22  1 2 

—  »  0-1*1+  X —  ^0^  plane  stress. 

A33  •'  A33 

(2.15) 
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•3. 


DISPLACEMENT  AND  STRESS  FIELDS  FOR  STRIPS 


The  two-dimensional  composite  medium  is  formed  of  two 
sets  of  periodically  arranged  strips  having  widths  2h]  and 
2h2  as  shown  in  Figure  1.  They  are  perfectly  bonded  along 
their  straight  boundaries,  and  contain  symmetrically  lo¬ 
cated  cracks  normal  to  the  interfaces,  of  length  2a  and  2b 
respectively.  The  load  is  applied  away  from  the  crack 
region,  such  that  the  crack  plane  is  a  plane  of  symmetry. 

Using  the  usual  superposition  technique,  the  solution 
of  the  actual  traction-free  crack  problem  may  be  obtained 
by  superposing  the  homogeneous  uncracked  strip  solution  to 
the  solution  of  a  cracked  strip  loaded  with  self-equilibra¬ 
ting  crack  surface  tractions  (see  Figure  2).  Since  we  are 
Interested  only  in  the  computation  of  stress  intensity 
factors  and  the  singular  behavior  of  the  stresses  around 
crack  ends,  we  will  consider  only  the  singular  part  of  the 
solution,  where  the  self-equilibrating  crack  tractions  are 
the  only  external  forces. 

First  we  will  find  solutions  to  (2.14)  satisfying 
certain  boundary  conditions  of  the  strips.  The  combination 
of  these  solutions  will  be  forced  to  satisfy  the  remaining 
boundary  and  continuity  conditions. 
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Assume: 


u^®^(x,y)  =  f  [  f(a,x)cos  ay  do 
Jq 

v^*^(x,y)  =  f  f  g(a,x3sin  oy  do 

JQ 

Substituting  (3.1)  into  (2.14)  we  obtain: 

S]  -  a"f  +  63  a  ^  =  0 

^  -  02  «"g  -  ^3  a  II  =  0 
The  solution  of  (3.2)  can  be  written  as: 

f(a,x)  =  A(o)e®‘“’'  +  B(o)e'*‘“’‘  +  C(o)e®"“’'  +  D(o)e 
g(o,x)  -  07[A(o)e®‘®''  -  B(o)e*®‘“’']  +  0g[c(o)e®^“'‘ 

-  D(a)e'®*“’‘] 

where  Sj  and  sz  are  the  roots  of 

s"  +  B^s*  +  Bg  =  0 


and 


B3  -  0,02  -  1 


B, 


“  0,s 


1  -  0,sf 


3^1 


0,s 


85  =  07  •  ^6  =  -^8^^ 
1®2 


1  -  B,s? 


3^2 


(3.1) 


(3.2) 


Szox 


(3.3) 


(3.4) 
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From  (3.4)  we  can  write 


and  $2  are  both  real  or  complex  conjugates. 


Assume: 

u^*’^(x,y)  =  f  [  h(o,y)sin  ax  da 
v^^’^Cx.y)  =  I  jo  A(o,y)cos  ax  da 
Substituting  (3.6)  into  (2.14)  we  have: 

0  -  63  a  ^  -  B^a^  h  =0 

h  0  ®  0  =0 

Solving  (3.7)  we  obtain: 

h(a,y)  «  E(a)e®‘®y''’^  +  F(a)e"® 

+  G(a)e^"®y/’^  +  H(a)e'®"“y^’^  , 

and 

t(a.y)  =  69(E(a)e®»®y'’^  -  F(a)e'®*“J'/’^J 


(3.5) 


(3.6) 


(3.7) 
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P,o[G(a)e**“y''^  -  H{a)e-®^“y'’^] 


where 


10  0, 


Bi/Bs  ^  S'! 


ii!S  +  ii- 

"2 


(3.8) 


A  superscript  *  will  be  used  for  the  material  constants  and 
unknown  functions  when  the  above  expressions  are  used  for 
the  second  strip. 

If  one  examines  the  roots  of  equation  (3.4),  he  will 
realize  that  there  are  two  types  of  orthotropic  materials. 
He  will  denote  the  material  as  type  I  when  s^  and  s^  are 
real,  and  as  type  II  when  they  are  complex  conjugates.  We 
will  assume,  in  our  analysis,  that  the  materials  of  both 
strips  are  of  type  I.  Similar  analysis  can  be  done  for 
the  remaining  combinations. 

3.3  The  Displacements 

for  each  strip,  we  can  write: 

u(x,y)  =*  u^*^(x,y)  +  u^**hx,y) 

v{x,y)  *  v^®^{x,y)  +  v^*’^(x,y)  . 

Noting  that; 

u(x.y)  =  -  u(-x,y)  and  v(x,y)  =  -  v(x,-y) 
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we  win  obtain: 

B{a)  =  -A'(a),  0(a)  =  -C(a),  F(a)  =  -E{a),  H{a)  =  -  G{a). 
For  material  type  I; 


Using  the  information  given  above,  and  keeping  in  mind  that 
u  and  V  vanish  when  y  goes  to  infinity,  for  y>0,  the  dis¬ 
placement  expressions  can  be  written  as  follows: 

u(x,y)  =  ^  |°*[A(a)sinh(w^ox)  +  C{o)sinh(w2ax)]cos  ay  da 


+  ^  f  E(a)e'l''‘  +  G{a]b->”3 


JoL 


sin  ax  da 


r 


v(x,y)  »  -  j^{B7A{a)cosh(wyax)  +  ggC(a)cosh(w3ax)]sin  ay  da 

-  f  |^|sign(W|)3gE(a)e’l*'*l°®'^^  +  sign(w3)G{a)B^(^'^*'^l®^'^'^ 


•cos  ax  da 


(3.9) 


3.4  The  stresses 


For  generalized  plane  stress  case: 

“  *n®x  *12^y 
®y  “  *12'^x  *  ^22^y 

■'xy  “  *33''^xy  •  (3.10) 
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Differentiating  (3.9)  and  using  (2.12)  and  (3.10), 
the  stress  expressions  can  be  written  as: 


2E,  -X 


(x,y)  =  r 


YiE(a)e 


-|wi|ay//g7  ^ 


Y2®(“) 


.-Iwslay/v^ 


ocos  ax  da  +  |[2Y2A(o)cosh(w^ax)  +  2Y4C(a) 


*cosh(w2ax)]acos  oy  da 


2Ey  ''y 


Y5E(a)e'*”‘l“y/'^+ YgG(a)  . 


(x.y)  =  f 

'**3l°‘y/'^iacos  ox  da  +  j^2Y7A(a)cosh(w^ox)  +  2YgC(a) 


cosh(w2ax)]acos  oy  da 


ZGjjy  Tjjy(x,y)  =  J  [2YgA(a)sinh(w^ox)  +  2y^o^^“^  • 

•  sinh(w2ax)]asin  ay  da  + 


(3.11) 

These  expressions  are  valid  also  for  the  plane  strain  case 
with  the  following  substitutions: 

'’yx  “  *12^*11  •  ''xy  “  *12/^22  ’  °  ^/*11  * 

(Ex*A)  =  I/A22  • 

The  elastic  material  constants  yj  are  defined  in  Appendix  A. 
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Yi2®(“) 


,-|w3|ay/y^ 


asin  ox  da 


Y„E(a)e-|^3|ay/^ 


4.  FORMULATION  OF  THE.  PROBLEM 

The  solution  of  the  problem  may  be  obtained  by  de¬ 
termining  the  unknown  functions  which  appear  in  the  dis¬ 
placement  and  stress  expressions,  under  the  following 
boundary  and  continuity  conditions: 

u^Ch^.y)  =  U2(-h2.y)  . 

v^(h^,y)  =  V2(-h2.y)(0  <  y  <  «>)  (4.1a,b) 

T^xyt^l>y)  *  ■^2xy^‘^2’^^^°  <  y<")  (4.2a.b) 

u,(0,y)  =  0  ,  T^xy(®»y)  =  0  (0  <  y  <  <»>)(4.3a,b) 
U2(0,y)  »  0  ,  T2j(y(0,y)  =  0  (0  <  y<«)(4.4a,b) 

^lxy^*l’^^”  ®  • 

T2xy(x2,0)  =  0  ,  !x2l<h2  (4.5a,b) 

Oly(xi  .0)  =•  -p^(x,)  ,  Ixy|<a  , 

v^(x^,0)  =  0  ,  a<iX||<h^  (4.6a,b) 

^2y^^2*^^  *’"P2^^2^  *  * 

V2(x2,0)  =  0  ,  b<Ix2l<h2  {4.7a,b) 

The  conditions  (4.3a,b)  and  (4.4a,b)  are  satisfied  iden¬ 
tically. 
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Using  (4.5a,b)  we  obtain: 

G(a)  =  -  E{a)  .  G*(a)  =  -  ^  E*(a) 


Ylz 


The  mixed  condition  (4.6)  gives: 


Tim  feCo) 


5  ''6  Y,j 


ocos  ax^da 


f  [2Y7A(a)cosh(w^  ox^)  +  2YgC(a)cosh(w2ax^ )]acos  ay  da 


2Ey^  ’  1*1 


(4.8a) 


and 


v^(x^,0)  =  -  I  Y^3|^E(a)cos  ax^da  =  0  a<|x^|<h^ 

(4.8b) 


Oefi ne , 


3v,(x, ,0) 

- g^q -  =  <(i(X|)  such  that  i}i(x^) 


=  0  for  1  x^  1>  a . 

(4.9) 


Differentiating  (4.8b)  with  respect  to  x^  and  taking  the 
inverse  transform,  we  obtain: 

Yi3aE(a)  =  |jj(()( x^ ) sin  ax^  dx^  .  (4.10) 

If  we  now  substitute  (4.10)  into  (4.8a)  and  evaluate 
some  of  the  integrals  in  tlosed  form  (see  Appendix  C)  we 
will  end  up  with  the  following  singular  integral  equation: 
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^14  I  dt  +  |”[2YyA(a)cosh(w^ax^ )  +  2YgC(a)cosh(w2a'«^  )]ada 


Pl^) 


-a<x^  <a 


(4.11) 


where  because  of  symmetry  4i(t)  =  -  (ti(-t), 


Similarly  defining, 

3v2(x2,0]  *  * 

— 331 -  =4i(x2)  such  thati)>(x2)  =  0  for  |x2]>b 

(4.12) 

and  using  the  mixed  condition  (4.7)  by  the  same  procedure 
we  obtain: 

w  <b*iti  ★  *  ★  ★  ★  ★ 

Yi4  j  '  dt  +  J  [2Y7A  (a)cosh(w^aX2)  +  2YgC  (o)cosh(w2aX2)  ]ada 


if(l-v*  V*  ) 


P2(X2)  -b<X2<b 


(4.13) 


The  next  step  is  to  determine  the  unknown  functions  A(a), 
C(a),  A*(a) ,  C*{a).  This  can  be  done  by  using  the  contin¬ 
uity  conditions  (4.1a,b)  and  (4.2a,b)  and  taking  the  in¬ 
verse  transforms.  Then  we  obtain  the  following  system  of 
linear  equations: 


2A(o)sinh(w^ah^ )  +  2C(a)sinh(w3ah^ j 
+  2A*(a)sinh(w*ah2)  +  2C*(a) si nh (Wjohg )  =  (a) 

2ByA(a)cosh(w^oh^ )  +  26gC(a)cosh(w2ah^ ) 

-  2B7  A*(a)cosh(w;|ah2)  -  2BgC*{a)cosh{w2ah2)  =  R2(“) 
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2Y3A(a)cosh(w^ah^ )  +  2Y^C(a)cosh (w^ah^ ) 

-  2X^T2A*{o)cosh(Wjah2)  -  2X^Y4C  {a)cosh(w2ah2)  =  R3(“) 

2YgA(a)sinh(w^ah^ )  +  2Y^QC(a)sinh(w2ah^ ) 

+  2YgX2A*(a)sinh(wjah2)  +  2YiQX2C*(a)sinh(w3ah2)  = R4{ct) 

(4.14) 


R|{a)  and  Xj  are  defined  in  Appendices  B  and  A  respectively. 
Solving  (4.14)  we  obtain: 

1  R-i  ( Ro (^) 

=  2cosh(w,ah,}  TT^  92'-®^  *’2^“^ 

R3(a)  R,{a) 

*  TT^  '"2(“)  *  TTST  "2^“^ 

R,{a)  RoCa) 

9,  (a)  +  h,{a) 


Ri(a) 


R4(a) 


“f(^  "»,(«)  +  ni(«) 


.  -if  -1  (®)  ® ) 

^  ^  2cosi;[w>2y  W  M®)  ^ 


RjCa) 


R4(«) 


*  TT^  "o<®^  *  -fwr 


Ri  (a) 


R,(ci) 


*  1  ■'ovu; 

C  (a)  *  - g(a)  +  -fr-\  h{a) 

2cosh(w3ah2)  ^^®' 


R3(a)  R4{a) 

*  "'^®^  ■"  Twr  "^®^i 


(4.15) 
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The  functions  f(a),  g(a)  fete.,  used  above  are  given  in 
Appendix  B. 

Substituting  (4.15}  into  (4.11)  and  (4.13)  we  obtain 
the  following  system  of  singular  integral  equations: 


?  r  ^  k,,(xpt)i(t}dt 


-  a  <  x^  <  a 


*  * 

-b 


-  b  <  x^  <  b 


(4.16a,b) 


where 


+  k^(x^.a)e’^V^^®'et/l'^3'  da 
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k2,(x2*t)  = 

14  °L 

+  kg(x2.a)e‘^V‘^“’^/l'"*l  da  (4.17) 

The  expressions  k.(j  =  1,8)  are  given  in  Appendix  B. 

J 

By  letting  h2-»*  (or  +  ®)  one  can  recover  the  spec¬ 
ial  case  studied  in  [5].  For  a  <  and  b<h2,  the  in¬ 
tegrands  of  kernels  k^j(i,j  =  1,2)  vanish  when  a-*-®  and  are 
bounded  for  all  values  of  a,  except  when  a  *  0. 

Around  a  =  0  the  asymptotic  behavior  of  the  inte¬ 
grands  I.,  of  the  kernels  k..  is  of  the  following  form: 

1  J  .  •  *  J 

^ij(a)  *  ^  ^  OH)  (i.j  =  1,2)  (4.18) 

where  the  c.-'s  are  known  constants.  In  order  to  obtain 
*  J 

a  solution,  one  should  show  that  the  singularity  due  to 
1/a  is  removable.  Consider  the  following  integral: 

I*  «(t)dt  |”l^^(x^,t,a)da  =  |  (j)(t)dt  (x^  ,t,a)da 

*  n  ^11  (x-j  ,t,a)da 
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where  e  is  a  positive  small  number.  Using  (4.18)  for  the 
first  part  of  the  integral,  we  obtain: 


I  (x^  ,t)({i(t)dt  =  I  4i(t)dt 
-a  -a 


r  ill 

Jo  “ 


da 


+  I  0(l)da  +  I  ^ (x^ ,t,o)da 


Making  use  of  the  single-valuedness  condition 

[  (jCtjdt  =  0,  the  unbounded  integral  [  da  drops  out, 

ia  Jo  “ 

leaving  only  bounded  integrals  which  can  be  evaluated  numer- 
* 

ically.  Using  <j)(t)dt  =  0  for  the  second  crack,  similarly 
ib 

one  can  show  that  the  singularity  due  to  1/a  cancels  in 
all  the  integrands  I... 

I  J 
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5. 


CASE  OF  BROKEN  LAMINATES 


This  is  the  case  when  one  of  the  cracks  touches  the 
interfaces  (i.e.,  a  =  h-j ,  or  b=h2).  The  integral  equa¬ 
tions  (4.16a,b)  are  still  valid  but  some  of  the  kernels 
kij  are  no  longer  bounded.  For  example  for  a =  hy  k^g* 
^21’  *^22  bounded  but  k^^  becomes  unbounded  as  and 
t  approach  the  ends  +  h^  simultaneously.  In  this  case 
the  integrand  I^^  of  k^^  diverges  as  o-»<». 

In  order  to  obtain  the  proper  singularity  at  the 
crack  tips  and  to  compute  k^^  numerically,  the  singular 
part,  should  be  evaluated  in  closed  form.  In  this 

case,  the  kernel  k^^  can  be  written  as: 


kii(xi,t)  =  ^lls^^l*^^  *  ^llf^*l’^^ 

where  k^^^  is  the  singular  part  and  k^^^  is  the  bounded 
part  of  k^  ^ . 

Following  the  procedure  described  in  [4]  k^j^^Cx^.t) 
is  obtained  as  follows: 

,  ,  ,  f  (hTt)vB7/lw,I  +  lwi|h, 

irkiis(Xi,t)  =  ^g5|^[(h^.t)/g7/|„jl  +  lwilh,]^-(w. 

((hi-t)vf7/Iwil  +  Iwsih, 

[(h^-t).fr/|w,|  +  |w3|h^]^-(w. 

^  ^  f  (hTt)v^/|w3|  +  |wi|h^ 

[(hTt).^/|W3|  +  |w3|h,]^-{Wi 


-25- 


+  X 


88 


((h-|-t)i^/|w3|  ■»  Iwalh^ _ 

[(h^-t)/eiy|w3|  +  iwalh^]^  -  (wjX^)^  ^ 


-  ,t  <  h-j 


(5.1) 


The  governing  singular  integral  equations  become: 

fh. 


1  1 


+  ’Tk^l5(xj,t) 


^  J  (t)dt  ~  2y^^£  ^  P-j  (x-j )  -  hi<  x-j<  h-j 


"1 

C 

rb 

lb  ''z 


<t>{t}dt+  [  ’ 

\l 


4»{t)dt 


J:  f*’  ^Z^dt  +  f^’k2,(x2,t)4i(t)dt  +  f**  k22(x2.t)4.*(t)dt  = 

L  9  k 


*  * 

»  V-  P2(X2)  -b<X2<b  (5.2a,b) 

^14^y 


Since  in  the  integral  equation  {5.2b)  the  only  singular 

term  is  -r-^r— ,  the  power  of  singularity  at  the  end  of  the 
t-x2 

internal  crack  in  the  second  layer  is  still  1/2.  But  in 
(5.2a)  we  have  further  singular  contribution  from  the  ker- 
ne1  k^^g  resulting  in  a  power  different  than  1/2.  To  find 
this  singularity  power  y.  we  will  again  use  the  procedure 
described  in  reference  [14].  throwing  all  the  bounded 
terms  to  the  right  hand  side,  the  singular  integral  equa¬ 
tion  (5.2a)  can  be  written  as  follows: 


^^11 s^^l 


<ti(t)dt  =  P(x,) 


-  hi<  Xi<  h^ 


(5.3) 
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where  P(x^)  is  a  bounded  function  for  all  values  of  . 

The  unknown  function  4i(t)  can  be  written  as  (see  [14]): 


(hf-t*)Y 


(5.4) 


where  F(t)  is  bounded  and  Holder-continuous  in  the  inter¬ 
val  |t|<h^,  and  0<Re(Y)<l- 

Define  the  sectionally  holomorphic  function: 


,  i  [^1  iLti  dt  =  1 


Fftie' 


Then , 


h^  (t-h^)^(t+h^)Y{t-z) 


F(.hi)e'"^  1  1 

- L -  — - - ! - ! -  + 

(2h^)^sinirY  (z+h^)^  (2h^)^siniTY  (z-h^)^ 

(5.5) 


where 

l>l'o(z)l<-; — ^-77-  »  Yo<Re(Y) 
Iz+h^To  ° 

C  and  Yq  SI'S  real  constants. 

Using  (5.5),  equation  (5.3)  takes  the  form: 


F{-hi)cotTrY 

(2h^)^(h^+x,)^ 


F(h^)cot7nr  Agg  |w^| 

(2h,)^(h,-x,)^  2 


F(h^) 

(w^  IT 

(2h,  )TsimiY 
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(h^+x^)^  (h^-Xi)^ 


^86  f”*! 
■7 


F(h,) 


( 


{2h^)^simrY 


1  1 


Iwill* 


w,iF 


/Bl 


(h,+x^)'F  (h^-x^)^ 


^87  I*',! 


F(h,) 


(2h^)^simrY 


Iwil jwal 


J - +  1 


(h^+X|)^  (h^-x^)"'^ 


+  - 


^88  l***^!  F(hi) 


{2h^)^simrY 


w 

^1 


1  +_L 


(h^+x^)"^ 


-  P,{x^) 


(5.6) 


where  because  of  symmetry  F(h.|)  =  -  F(-h^). 

Multiplying  both  sides  of  (5.6)  by  (h^+x^)^  and  let¬ 
ting  x-j  =  .-  h-j  we  obtain  the  following  characteristic 
equation: 


-  2  cosiTf  +  Agg 


KI 

1 

|WiI 

„  4.  1  _  _ 

1 

/b7 

KIIW3I 

Y 

L  J 

'‘S? 


IM 

1 

|W-| 

.a.  -1  1  3  1 

1 

Kih,i 

«3 

Y 

.  ^  . 

=  0 


(5.7) 


where  are  elastic  constants  defined  in  Appendix  A. 
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This  Is  the  same  equation  found  in  [7].  Choosing  the  or¬ 
thotropic  elastic  constants  close  to  isotropic  constants 
numerically  we  find  the  same  singularity  power  computed  in 
[6]  and  [9].  The  characteristic  equation  (5.7)  can  be 
solved  numerically  to  find  y.  For  practical  orthotropic 
materials  equation  (5.7)  has  only  one  root  between  0  and  1. 
To  establish  the  dependence  of  y  on  the  material  constants 
more  accurately,  a  separate  study  of  equation  (5.7)  is 
needed. 
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6.  CASE  ^  A  CRACK  CROSSING  THE  INTERFACE 

To  formulate  this  problem  we  will  start  by  using  the 
crack  configuration  shown  in  Figure  3.  In  this  case  we 
have  again  an  internal  crack  in  the  first  layer,  but  two 
symmetrically  located  cracks  in  the  second  layer.  Using 
the  symmetry  property  of  i{i*{t),  we  can  write: 


I  k^2{Xi>t)'i'*(t)dt  =  ■  '‘i2(>‘i»-t)]<t'*(t)dt(i  =  1,2) 

•b 


and 


r  ^ •  f' 

lo  2  '• 


f*’ 

1 

^  1 

Jo 

t-X2 

(J>*(t)dt 


(6.1) 


Therefore  we  can  write  the  governing  singular  integral 
equations,  by  simply  changing  the  limits  of  the  integrals 
from  (0,b)  to  (c,d)  in  equations  {4.16a,b).  Thus, we 
obtain: 


ir^'it"rhi(’‘i>t)^(t)dt 


’  ■  Pi^"1^  -a<x^<a 


1 

i  c 


1  +  1 


t-X2  t+Xg 


14‘-y 
*  f® 

<(i  (t)dt  +  J  k2^(x2,t)ii)(t)dt 
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+  J  ° - -..y^y  p^(x^) 

c  Zy^^Ey 

c<X2<d  ,{6.2a,b) 

By  letting  a  =  h^  and  d  =  h2  we  obtain  the  case  of  a  crack 
crossing  the  interface.  As  in  the  previous  case  for  a  =  h^ 
and  d  =  h2  all  the  kernels  k^.j  become  unbounded  when  x^,t 
approach  the  ends  +  h-j  and  X2,t  approach  the  end  h2  simul¬ 
taneously.  Therefore  to  study  the  singular  behavior  at  the 
interface  and  to  make  the  kernels  numerically  integrable, 
the  singular  parts  of  the  kernels  k..  must  be  separated. 

1  J 

The  kernels  can  be  written  as: 

k^j(x^,t)  =  k^jg(x^,t)  +  k^jf(x^..t) 

where  the  singular  and  k^j^rCx^.t)  is  the 

bounded  part.  Following  the  same  procedure  used  in  the 
previous  section  the  expressions  of  k^^^  are  found  as 
fol lows: 

ffklls(Xl,t)  =•  Xgg 


E(h^-t)i/B7/|w3|  +  |wi|h^  "1 
[(h,-tyBF/|w,|  +  lwi|hi]^-{wiX^)^J 

t(h^-t)yB7/|w3|  +  |w3|h^  I 
[(h,-t)y8i/|w,|+|w3|h^]^-{w3X^)^) 


f  (h^-t)/^/Iwi|  +  |w3|h^ 

®®|^C{h^-t)K^/|wi}+|w3lhi]^  -  (wjx 


{ruvt 


-t)yB7/|wii 


l  +  |wi|h^  ^ 

iJ^-(vi)7 
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rt,2,(x,.t) .  A, 3! — *  i"'i''i _ ) 

(C(h2-t)»^/|wT|  +  lwi|h^]^  -  (w^x^)2j 

(h2-t)i/^/|Wi|  +  Iwalh^ 


* 


7') 


95 


+  X 


*  X 


C(h2-t)i/^/|wi  I  +  Iwalhj]^  -  (w^x 

f  ■'■  l'*'i  1^1  ) 

ic{h2-t)i^/|wtl  +  Iwajh^]^  -  {w^x^)^] 

j  (Vt),/^/|wth  |w3|h^  \ 

'\[{h2-t)i^/lw3l  +  Iwalh^]^  -  (w3X^)2j 

r  {h'j”t)ii^/|wi  I  +  |wi  I  hg 
I  +  |wT|h2]^  -  (w’^X2)^ 

I  (hi-t)i/^/|wi  I  +  |wt|h2 
(  [(hj-t)./f7/iwil  +  Iwtjhg]^  -  (wjXg)^ 

+  X  (  (hi-t)/B?/lw3|  ■»  Iwtihg  I 

(  [(h|-t)»^/|w3|  +  |wt|h2]*  -  (w^X2)*J 
(h'j-t)/57/lw3|  +  |w3|h2 
■t)i^/|w3|  +  Iwjjhg]^  -  ((11(3X2)^ 

(h2-t)v^/|wi|  +  |wt|h2  \ 
C(h2-t)i^/|wt|  +  (wt(h2]^  -  (w*X2)^J 
(h2-t)i'^ /|wt|  +  |wt|h2  I 
[(h2-t)i/^/Iwt|  +  Iwtlh^]^-  (w3X2)^j 
{h2~t)*'^/|w3 1+  |wi  I  ^2 


irk22s(x2’*)  =  X^ggj 


*  ^110 1 

*  ^111  I 

*  "172  I 


[(h2-t).^/lw3!  +  Iwilhg]*  -  (w 
(h2-t)''^/|w3 1  +  |wt|  ^2 

[(h2-t)i/^/|w3|  +  Iw3|h2]^  -  (W3X2) 
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Separating  the  singular  parts,  equations  {6.2a,b)  take  the 
form: 


f”! 

+  j  [k^^(x^,t)  -  j(x^  ,t|()){t)dt 
ihi 

+  |^^k*25(x^,t)<).*(t)dt 


*  I  ^C'‘*2^’‘l ■  k*2s(^l  >^)3't’*(t)dt 


— p^  (x^ )  -  h^<  x^<  h, 


2y 


and 


i  f' 

If  Jr 


14‘-y 


T'iq  *  tnq  * 


{t)dt 


^  J  ”  1^2 1  s  ^ ^2  ’ ^  J  J ^ ^  J 

-h, 

fh-i 


t^2  * 

••c  ^*‘22^*2’^^  -  k22s(x2’^^ 


*  * 

(!.'»  ''  ) 

-  .  P2(x2)  c<X2<h2 


(6.3a, b) 
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where 


k*2(x^.t)  =■  k.gCx-.t)  -  k.2(x..-t)  (1  =  1.2)  . 

To  find  the  proper  singularity  power  B  at  the  inter¬ 
face,  we  will  first  throw  all  the  bounded  terms  to  the 
right  hand  sides  of  the  equations  as  it  has  been  done  in 
the  previous  case.  Then,  we  obtain  the  following  system 
of  equations: 


Ti  \  [t-V  *  irk^^j(x^.t)  i))(t)dt 

-“i  L  . '  J 


1  (**2  r  1 

7Jc  [t-: 


-h^<x,<h^ 


^2  *  ^^^^2  *  ’^*^225^^2’^^ 


fh, 

-h, 

(6.4a.b) 


4'*(t)dt 


1 


c<  X2<  h2 


where  Q^Cx-j)  and  Q2{x2)  are  bounded  functions  of  x-j  and  Xg. 
Considering  the  behavior  of  <{>(t)  and  ()i*(t)  at  the  end 
points,  we  can  write: 


<(>{t)  = 


F(t) 


<C*(t)  = 


F  (t) 


(h2-t)e(t-c)« 


(6. 5a,b) 


(h^-t^)'’ 

Define  the  following  sectionally  holomorphic  functions 
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■s 


From  [14]  and  using  (6.5a,b)  and  (6.6a,b)  we  have: 


(2h^)®simr6  (z+h^)®  (2h|)^siniTB  (z-h^ 


6  +  'f'o(z) 


,-(,1  .  -£!is]^ — 7^-5**;  w 

(h2-c)^sirro5  (z-c)*  (h2-c)°simre  (z-h2)® 

(6.7a.b) 

where  <l'g(z)  and  <I<q{z)  are  bounded  functions  which  around 
ends  behave  as  follows: 

C_  . 


Iz+h^l^o 


B„<Re(B) 


and 


|2-h2l 


Iz-cl^o 


^  .  Bo<Re(B) 


S„<  Re(6) 


*'0*  '^o’  ^o’  ^o’  ®o’  *0  **’®  constants. 

Using  (6.7a,b}  and  following  the  procedure  used  in 
section  5,  equations  (6.4a,b)  reduce  to: 


cotird  =  0 


(6.8) 
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and 


F  (h-) 

+ - 1 - 

2(h2-c)  sinirB 
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1*1  iu*l  1 8  112  ^  wt^  1® 


(6.9a.b) 


Equation  (6.8)  gives  the  expected  6  =  1/2  singularity  power 


at  the  crack  tip. 


(6.9a,b)  is  a  system  of  homogeneous  linear  equations 
for  F(h^)  and  F*(h2).  Since  F(h^)  /  0,  F*(h2)  f  0, 

6  f  0,1  to  solve  the  system  one  should  equate  the  determ¬ 
inant  of  coefficients  to  zero.  Thus, 


A(8)  »  4cos*ir6  -  Zcosirg 


fx  Ki-J- 

Vl09  j-r  (*  1 


^  ('•i/  (’•33)^ 


.1  i-ii  1 

^  ^  (r„)  ® 


+  1  A  1  -- 


+  X  _ 

^  (r33)  ^ 
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+  X  1  +  X  1 

*115  (  *  ,6  ^  *116  ,  *  .6 

1'^13'^iv 


X  1  +1  1  ,  ,  1 

*117  .  *  .6  *118  ,  *  ^  *119  ,  *  .8 

iriaPiaJ  1*^33*^111 


X  1  +  X  1 

*120  *  jB  *121  ,  *  jB 

1*^33*^131  1*^33*^331 


(6.10) 


where 


hiihsi 

‘11 

a  .  f 

’  *■33 

* 

w** 

.  1 

.  _*  .  Iwtllwtl 

if 

11 

*  -  » 

”  ■ 

’  '■33 

From  the  characteristic  equation  (6.10)  we  can  determine 
the  singularity  power  B.  -Choosing  the  orthotropic  elastic 
constants  close  to  isotropic  elastic  constants,  we  recover 
the  singularity  power  found  in  [8]  and  [9].  Equation 
(6.10)  does  not  always  have  a  root  between  0  and  1.  For 
some  orthotropic  material  combinations,  there  is  no  power 
singularity  at  the  interface.  In  this  case  one  should  in¬ 
vestigate  the  possibility  of  pure  imaginary  or  complex 
roots.  Numerical  computation  shows  that  there  are  no  pure 
imaginary  roots  or  complex  roots  for  which  the  real  part 


Is  between  0  and  1.  On  the  other  hand  F(h^)  and  F*{h2) 
are  related  through  (6.9a)  or  (6.9b>.  This  is  a  condition 
to  be  used  while  obtaining  the  solution.  The  absence  of 
power  singularity  for  some  orfhotropic  material  combina¬ 
tions  may  be  very  important  from  the  view  point  of  design 
applications.  Therefore  we  will  study  in  some  detail  the 
behavior  of  the  crack  surface  displacement  derivatives  and 
the  stresses  at  the  interface. 

Let's  first  investigate  the  possibility  of  a  weaker 
I.e.,  logarithmic  singularity  in  the  crack  surface  dis¬ 
placement  derivatives  at  the  interface.  Suppose  that  the 
power  singularity  0  at  the  interface  is  zero.  Define: 


The  behavior  of  ii/{z)  around  z  =  +  h^ ,  and  of  ii>*{z)  around 
z  =  h2  can  be  expressed  as: 

i|»*(z)  =•  log(z-hi)  +  l^oiCz)  near  z  =  h^ 

'f'(-hi) 

>1'  (z)  =  -  — -  log(z+h^)  +  near  z  =  -h^ 

*  ^  ( h? ) 

<|i  (z)  »  - log(z-h2)  +  'f>o3(z)  hear  z  =  +  h2 

(6.11a,b,c} 

where  ((|q^{z),  <(1^2(2),  '{'03(2}  are  bounded  functions. 
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Using  (6.11a, b,c},  near  = h^  and  Xg  =  hg  equations 
(6.4a,b)  take  the  following  form: 


(6.12a,b) 


where  F^(xi)  and  F2(x2)  are  bounded  functions. 

In  order  that  equations  (6.12  a,b)  be  bounded  for  Xj=h,,  x^s 

h2,  the  coefficients  of  the  logarithmic  terms  must  be  zero. 
Thus: 


4{h^)  |w^|  |wj  |wj  |wj 

^  2^  '  2^  ’  ^  2/b7 


|w*|  |w*|  |w*|  |w*| 

'  2^  ’  2^  '  2v^  '  2/^ 


and 


(6.13a,b)  is  a  system  of  linear  equations  for  iji(h^)  and 
<fr*(h2).  Since  and  <}i*(h2)  are  different  than  zero, 

in  order  to  have  a  solution  the  determinant  of  coefficients, 
A,  must  be  zero.  Numerical  computation  shows  that  A  =  0 
and  either  from  equation  {6.13a)  or  (6.13b)  we  have: 

<|i  (h,) 

«  '  =  -  1 

4  (h2) 

•ff  tff 

and  using  the  symmetry  condition  ip  (h2)  =  -  p  (-hg).  we 
obtain; 


♦  (h,) 

/(-hg) 


(6.14) 


Relation  (6.14)  shows  that  the  surface  displacement  deriv¬ 
ative  is  continuous  at  the  interface.  This  is  an  important 
result  which  makes  the  solution  of  the  singular  integral 
equations  easier. 


To  study  the  behavior  of  the  stresses,  let's  first 
write  their  expressions  at  the  interface.  By  making  use  of 


(3.11)  and  (4.15)  we  obtain: 

a{y)  s  {h|,y)  =  |jic^(n.y)$(n)dfi +  |^^K2(n.y)i{i*(n)dn 

T(y)  »  ^  ■'lxy^*’l’^^  '  |**^K3(n,y)$(n)dn  +  |  ^K^(n,y)il'*(Ti)dn 

(6.15a,b) 


where  K.(n,y)  (j  =  1,4)  are  given  in  Appendix  B. 


The  kernels  Kj  become  unbounded  as  y-»0*  and  n -*•  i  h^ 
or  n  ♦  h2  respectively.  When  equation  (6.10)  has  a  root, 
i.e.  when  the  functions  <))(t)  and  .4*(t)  are  singular,  the 
stresses  have  the  same  singularity  power  as  we  will  see  in 
the  derivation  of  the  stress  intensity  factors.  But  it  is 
necessary  to  know  whether  the  stresses  have  a' logari thmic 
singularity  when  the  crack  surface  displacement  derivatives 
are  bounded.  To  do  that  let's  first  separate  the  singular 
parts  of  kernels  K..  Following  the  usual  procedure,  we 
obtain: 


(h^ -n) 


Y]  (h^+n) 

W^V^  2  ^1  2 

-^+(hi+n)  -^+(hj-n) 


Yii  Tn 


+  Y,  — -Y-I- 
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IC3s(n.y)  = 


K45(n.y)  = 


1  (h,-n)./5i7|w  I 

■  /—  \  A-'-  \  \  ’  1 


*  2^I3»80  ‘''3*8'*^'''='’  [(l.,-n)«|..|]‘*y' 


1  ,  ^  (h,-n)*^/|wj 

ZY13X30  3  83  4  84  [{h  .n)^/lw3l]*+y^ 


1 

*  3  ^89%^90^  r/.  _3  r-r,i.  *.-.2  ,  ..2 


2^13^80 


^  (Yo^fli'^ii^qo)  _  *'  2  2 

2^13^80  [(h2-n)>^*/|wj]  +y 


[(h2-n)^/|w*|]  +y‘ 
(h2-ii)./^/lw*| 


Yl1 


|wi|y/i/^  |wiIy/»^ 


■2  2 
Wiy 


-r  •••  -z-z- 


L^Mh,«)  =g-*(l.,-nl 


+  — , 


Iwaly/ZST  Iwjly/*^ 


rci" 


w*y*’  i  w^y*  2 

-^+(h,+n)  -f7-+(h,-n) 


*  27j^A^  ^■'^9^81^10^82^ 


y 


C(hi-n)/8F/|w,13"  +  y' 


*  2^  <Vs3*3,0>84' 


2  2 
+  y 


^^13^80 


^■'^9^89'^^10^90^ 


[(h2-n)/^/|wt|]^  +  y* 


+  *  tyq^qi‘^in^q?5 


^.3^  V.,^,  .,0-33- 


(6.16a,b,c,d) 
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Keeping  only  the  singular  terras,  equations  (6.15a,b)  can 
be  written  as  follows: 

o{y)  »  I  ^K.|j{n,y)<t'(n)dn  +  |  ^K2j(n.y)it>*{n)dri  +  A{y) 

T(y)  »  I  ^K2j(n,y)<ti(n)dn  +  |  ^K2j{ri,y)i)>*{n)dn  +  B(y)  (6.17a,b) 

where  A(y)  and  B(y)  are  bounded  functions. 

Define: 

^(z)  =  I  1  dt  and  i|)*(z)  =  I  ^  dt 
-h, 

Considering  their  behavior  around  ends,  we  can  write; 
i|)(2)  »  <|i(h^)log{z-h^)  -  i{i(-h^)log(z+h^)  +  i(i^(z) 


and 


i|)*(z)  =  <(i*{h2)log(z-h2)  +  i(i*(z) 


(6.18a,b) 


where  'l'jj(z)  artd  ’)'q(z)  are  bounded  functions.  Making  use 
of  (6.18a,b)  and  (6.14),  equations  (6.17a,b)  become: 


a(y)  =  log  y  4i(h^) 

1  Iwj 


^''^13^80  ^ 

1  Iwtl 

^■*^13^80 


J_  (y  -y  Ill)  .  —I _  (y  X  +Y  X  ) 

'^13  ^2  y,2^  2Y13A80  ^  ^Vsi  ^4^82) 

I  w*  1 

■  (t^3^83%^84^  *  r*T  TT"  ^^3^89%^90^ 


^^3^9l'*”’^4^92^ 


+  C(y) 


-44- 


isti 


and 

T(y)  =•  D(y)  (6.19a,b) 

where  C(y)  and  D(y)  are  bounded  functions.  Equation 
(6.19b)  Indicates , that  the  shear  stress  at  the  interface 
Is  bounded.  Since  C(y)  is  bounded  and  (fr(h^)  f  0,  If  a(y) 

Is  bounded  the  coefficient  of  logy  In  equation  (6.19a) 
should  be  zero.  Numerical  computation  shows  that  the 
above  mentioned  coefficient  Is  Identically  zero.  Therefore 
the  normal  stress  o^  at  the  interface  is  also  bounded. 

These  are  important  results  for  orthotropic  materials  and 
may  have  practical  implications  in  designing  with  composite 
materials. 


THE  SOLUTION  AND  THE  RESULTS 


Since  we  have  mainly  three  different  problems,  the 
solution  will  be  discussed  in  three  sections. 

7.1  Case  of  Internal  Cracks 


In  this  case  we  have  to  solve  the  system  of  singular 
equations  (4.16a,b).  Defining  x.|  =  aKj ,  t  =  ar  for 
-  a<  ,  t  <a  and  X2  =  aK2.  t  =  bt  for  -b<X2,  t<  b  after 
normalization,  equations  {4.16a,b)  take  the  form: 


i’ll  ®  |^k°,(lC,,T)<fr^(T)dT  +  b  |^k°2{<|.T)<>o(T)dT 


1 

irj,  T-K- 


dT  +  a 


:i 


,  fl  *  ■ 

k2i(K2.T)'tiQ(T)dT  +  b  k22(K2>T)i!>(j(T)dT 

-1 

(7.1a,b) 


ZYl4Ey 


where  the  index  "0"  denotes  the  normalized  quantities.  To 
get  the  complete  solution  we  need  also  the  single-val ued- 
ness  conditions: 


(7.2a,b) 


Since  $o{t)  and  have  a  power  singularity  -1/2  at 


the  ends,  the  solution  will  be  sought  in  the  form: 

F„(t)  *  F*(t) 

=*  °  and  =  — - - 


/l^ 


where  Fj,{t)  and  F^^Cx)  are  Holder  continuous  in  the  interval 

-  1  <  T  <  1 . 


Using  the  method  described  in  [11]  we  obtain: 

N  ""  N 

•  -  "  SyJ/ 

N  N  **  ” 

*  * 

”-N - ?^P2(<i)  i  =  1....N-1 

2Yi4Ey 

N  '  N 

.1^  ^  “  °  i,  IT  =  °  (7.3a,b.c,< 


where 

Tj  »  cos  ^  (2j-l)  j  =  1,...N 
1C .  =  cos  ^  i  =  1 ,  — N-1 

The  2N  unknowns  F|j(Tj)  and  F^Cxj)  can  be  found  by  solving 
equations  ( 7. 3a ,b,c ,d) .  In  this  problem  we  are  mostly  in¬ 
terested  in  the  computation  of  the  stress  intensity  factors. 
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The  stress  intensity  factors  may  be  expressed  in  terms  of 
the  density  functions  Fq(t)  and  Fjj(t)  as  follows: 


For  a < h, :  k,  =  lim  /2(x, -a)  o,„(x,,0) 

'  “  xi^a  '  '3  ' 

and  for  b<h,:  ku  =  Tim  i/2{x,-b)  o,  (x,,0)  (7.4a,b) 

4  0  X2^b  ^  ‘y  ‘ 

Making  use  of  equations  (4.16a,b)  and  definitions  (7.4a, b), 
after  lengthy  algebra  (see  Appendix  D)  we  obtain: 

2yi 

'‘a  “  ■  (1-v  V — T  ^0^’^ 

'  xy  yx' 


n-Vyx^ 


(7.5a,b) 


The  computation  is  done  for  generalized  plane  stress  case 
only.  Results  can  easily  be  obtained  for  plane  strain  case 
by  redefining  the  elastic  material  constants.  In  the  per¬ 
turbation  problem  considered  p^  and  are  constant.  As¬ 
suming  that  there  is  no  constraint  in  x-  direction,  p^  and 
Pg  satisfy  the  following  condition: 

P2  ■  Ej 

where  E^  and  E*  are  the  Young's  moduli  in  y  direction. 

Two  material  combinations  are  formed  among  the  follow¬ 
ing  three  materials. 
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1.  Ejj  =  55. let  10^  N/m^  .  =  17.0.8X  10®  N/m^  , 

Sjjy  =  4.82  X  10®  N/in^  .  =  0.036 

(E^  =  8  X  10®  psi)  .  (E^  =  24.75  x  10®  psi)  , 

=  0.7  X  10®  psi)  .  (v^y  =  0.036) 

2.  Ejj  =134.4  X  10®  N/m^  .  =  31  .OPt  10®  N/m^  . 

Gj^  =24,12x  10®  N//  ,  v^y  =  0.650 

(Ejj  *  19.5  x  10®  psi)  ,  (Ey  =  4.5  x  10®  psi)  . 

(Gjjy  =  3.5  i  10®  psi)  .  (v^y  =  0.650) 

3.  Ejj  =154.7  X  10®  N/m^  .  Ey  =155.8  x  10®  N/m^  , 

Sj,y  =59.65X  10®  N/ra^  ,  v^^y  =  0.300 

(E^  =  22.447  X  10®  psi)  ,  (Ey  =  22.6  x  10®  psi)  . 

(Gjjy  =  8.655  X  10®  psi)  .  (v^y  =  0.300) 

As  it  is  seen  from  the  values  given  above  the  first  two 
materials  are  orthotropic,  while  the  third  is  isotropic. 
The  following  pairs  of  materials  are  used: 

Combination  I:  The  first  layer  is  of  material  1, 
the  second  of  material  2. 
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Combination  II:  The  first  layer  is  of  material 
3.  the  second  of  material  2. 

Choosing  the  same  materials  and  letting  a,  b,  h^  or  h^  go 
to  proper  limits  we  recover  all  the  special  cases  done  in 
[5],  [6],  and  [10]. 

Figures  4-12  show  some  of  the  calculated  results.  In 
Figures  4  and  5  the  stress  intensity  factors  k,  are  plotted 

d 

versus  h2/h^  for  b  =  0  (there  is  no  crack  in  the  second  ma¬ 
terial)  and  for  the  two  material  combinations.  For  h2 = 0, 
we  recover  the  solution  of  col  inear  cracks  imbedded  in  a 
homogeneous  material  (see  [10]).  It  is  important  to  note 
that  in  the  colinear  crack  problem  the  material  doesn't 
have  to  be  isotropic.  As  h2  +  “>,  k^  reaches  an  asymptotic 
value  which  can  be  found  in  [5].  For  a  fixed  h2/h^  ratio, 
increases  as  a/h^  becomes  larger. 

Figures  6  and  7  show  the  stress  intensity  factors  k|j 
for  the  case  a  =  0.  In  this  case  also,  for  h^  =  0,  we  ob¬ 
tain  the  solution  of  colinear  cracks.  There  is  a  critical 
value  of  (h^/h2)  for  which  the  stress  intensity  factor 
starts  to  decrease  as  the  ratio  b/hg  increases.  For  the 
examples  done  this  critical  ratio  is  between  0  and  0.5. 

For  h^+<»  the  stress  intensity  factor  kj^  reaches  an  asymp¬ 
totic  value  which  also  can  be  found  in  [5]. 
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The  stress  intensity  factors  and  when  both 
layers  contain  cracks,  are  given  in  Figures  8-11.  kj^-^O 
as  (b/hg)-*-!,  since  the  power  singularity  y  is  less  than 
0.5  when  the  crack  in  the  second  material  touches  the  in¬ 
terface.  Another  interesting  result  is  obtained  from  the, 
comparison  of  isotropic  and  orthotropic  materials.  As  it 
is  seen  in  Figure  12,  for  the  same  and  E*  the  stress  in 
tensity  factor  k,  for  orthotropic  materials  can  be  larger 

a 

or  smaller  than  the  stress  intensity  factor  k,  for  iso- 

d 

tropic  materials  depending  on  the  other  elastic  constants. 
One  can  significantly  reduce  k  by  a  convenient  choice  of 

a 

the  elastic  constants.  The  materials  used  in  the  compari¬ 
son  are  given  in  Table  1.  The  dependence  of  k^  on  the 

* 

materials  constants  is  given  in  Table  2.  G  and  G  are 

xy  xy 

it 

the  most  important  constants,  while  keeping  E^  and  E^  con¬ 
stant.  To  reduce  the  stress  intensity  factor  k^,  it  is 

.  'ft  Hr  ^ 

sufficient  to  increase  E  ,  G  ,  v  or  decrease  E  ,  G  , 

A  *y  ^y  ^ 
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7*2  Case  of  Broken  Laminates 

The  solution  will  be  obtained  by  solving  equations 
(5.2a,b),  with  the  single-valuedness  conditions, 

f  ^(>(t)dt  =  0  and  f  $*(t)dt  =  ff  {7.6a,b) 

■hi 

Defining  again, 

t  -  h,T^  .  X,  =  h,<^  for  -h,  <x,,t<h, 


1  '  "1 


t  »  br^  ,  *2  “  *”^2  -b<X2,t<b 

the  normalized  form  of  equations  {5.2a,b)  and(7.6a,b)  can 
be  written  as  follows: 

+  b  {^k°2('CpT2)<fr*(T2)dT2  =  - 

1  n  n  * 

¥  ¥^^*^2  *  ^]^'‘21^*'2’''l^'‘’o^^l^‘*^1  *  bJ^k22^*'2’^2^'*’o^'^2^'*''2° 


*  ♦ 

- J  I  -  P2('<2^  »  -1<|C2<1 


rl  * 

]  ~  ®  ]  ‘*’o^^2^*^^2  ■  (7.7a,b,c,d) 


To  obtain  the  solution,  we  will  use  the  numerical  method 
described  in  [11].  Hence,  we  obtain: 


*  ^  “j^^lj^'‘llf(’'li*'^lj^''o(^lj^ 

N 

*  **  N"  j  =  l''l2^'^li’’^2j^''o^^2j^ 


^^■'’xy'^yx^  0/^  ,  T  -  1  Ml 

IyTJe/  Pl(  li)  ’  ’  - 


1  7  f— 1—  . 

^  j  =  l  ^2j-'^2i 


irbk22{<2i*^2j^  To^'^2j^ 


*  *'j^‘^lj^''21^'^2i’‘^lj^*'o^'^lj^ 


» - pO(K^.)  ,  i  =  1....Nt1 

N  N 


(7.8a,b,c,d) 


where 


*«'  '*  ■  (T:^  ’ 


PN^‘^’'^^Tij)  =  0  ,  j  =  1,...N 

PNiJ'^^^’^'A^li)  =  0  ,  i  =  1,...N-1 
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‘‘b  “  xy4)  O2y(x2.0) 

After  some  calculation  shown  in  Appendix  D,  we  have: 


^  (1-v*  V*  }  ’  (7.9a.b) 

'  xy  yx' 

The  results  for  the  case  of  broken  laminates  are  shown  in 
Figures  13-16.  Again  the  same  material  combinations  are 
used.  In  Figure  13  the  stress  intensity  factor  k,  is 

a 

plotted  versus  the  ratio  h2/h^  for  b  =  0.  When  h2-*-<«>  k^ 

has  an  asymptotic  value  which  can  be  recovered  in  [7]. 

Figure  14  shows  the  variation  of  kj^  with  h^/h2,  for  the 

case  a  =  0.  For  b  /  0,  the  variation  of  k,  and  k^  with 

a  b 

b/h2  are  given  in  Figures  15  and  16. 

7.3  Case  of  a  Crack  Crossing  the  Interface 

In  this  case  the  governing  singular  integral  equa¬ 
tions  are  (6.3a,b).  As  it  was  pointed  out  in  Section  6, 
the  characteristic  equation  (6.10)  does  not  always  give  a 
singularity  power  at  the  bimaterial  interface.  Therefore, 
the  numerical  solution  needs  care,  and  we  should  solve 
equations  (6.3a,b)  considering  the  singular  and  non-, 
singular  cases  at  the  bimaterial  interfaces. 

7,3.1  Singular  Behavior  at  the  Interface 

For  the  material  combination  II  (isotropic-ortho¬ 
tropic)  equation  (6.10)  has  a  root  between  0  and  1.  Using 
Newton-Raphson  method  to  solve  equation  (6.10),  we  have: 

B  =  0.04248  . 
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We  win  make  use  of  the  following  definitions  to  normalize 
equations  (6.3a,b),  the  singl e-val uedness  condition,  and 
the  relation  (6.9a): 


t  »  h^T^  ,  =  h^K^  for  -h^<t,  X|<h^ 


h^-c  h2+c  hp-c  h-+c 

t  =  t2  +  >  X2  =  <2  ~2~ 


Then,  we  obtain: 


1  T,-tc,  ‘*'^1  *  *’1  1  '‘ns^''i’^i^‘*’o^''P‘*^i  *  N  j  '‘iif^''r''i^'*’o^^i^‘*^i 

-1  ’  ’  -1 


A  n**C  fl  -it 

*  ~2~  J  *‘^2$ ]^'‘l2f^’'l’'^2^'^o^^2^*2 


h2-c  fl 

^2-^+-L‘ 


-0. 


<])  -1<X^<1 


T  [  *  TTTTlTv^  '^'o^^2^*^'^2 


1  ^2-''2  ■x,+K,+  ‘a^ 

I  2  2  h2-C 


**1  [^'‘21s^''2’'^l^‘*’o^'^1^‘*^l  *  *hf^''21f^''2’''l^'^o^‘'l^'^’^l 


f'l  .^Q  if  n«*c  fl  if 


ho-c  fl 


*  *  . 


-1<K2<1 


2v  E 

hi  =  0  h^  ^  ^ 


(7.10a,b,c,d) 
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wHere 


Using  the  numerical  method  given  in  [11]  equations  (7.10 
a,b,c,d)  further  reduce  to: 


i  .1  T,  -K, .  *  ’^*’l''lls^''li’''lj^  ’^'’l''?lf^'^li’'^lj^  ''o^'^lj^ 

I J  1 1 

J!  *  hp-c  *0  ,  ,  *0  ,  .  *,  V 

*  ^  “T” '^12s^''li’^2j^  ~F" ’‘iZf^'^li’^Ej^ 


i  =  l....N-l 
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CiipuiUi.gliUji 


wgaijpgiiipp 


In*  1  1  h2*^  *0 

^  L  ■'2j'^'"2i^  hg-c 


h2-c  *0 


+  ir^k22f(.C2^.T2j)jFo{T2j)  *  Wj(^lj)['’l'‘l2s(''2r-^lj) 

(1  V*  V*  ) 

+  °  r  P2(''2i^  ’  °  1....N-1 

^Yi4ty 


*'1  =  ° 


h,  ®  a^ 

F*(l)  =■  -  ;? 


(7.na,b,c,d) 


where 


*  F*(t) 

A  (-p)  3 -  ^  (x)  =  — - r: 

°  (i-i‘)6  »  (i-T)e(UT)'< 


3  =  1,....N 


p^^^0-6,1-B)(^^i)  =  0  1  =  1,...N-1 

p^(-6,-*s)(^^^)  ,  0  j  =  1,...N 

P^_^0-B.>s)(^2,.)  =  0  i  =  l....N-l 


Wj(T^j)  and  Wj(T2j)  are  the  corresponding  weights  of 
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p  (-6.-6)(  )  p 


'N  —  -N 

Solving  the  2N  x  2N  system  of  linear  equations,  we 
obtain  the  2N  unknowns  Pg(fij)  snd  F*(T2j).  The  stress  In¬ 
tensity  factors  can  be  defined  as  follows:  -  , 

kj,  -  11m  ✓2(c-X2)  O2y(x2,0) 

and  at  the  bimaterial  interfaces 


■  ,’Jo’  ' 

‘xy  ■  ,Vo"  'lx,"’!-'' 


(7.12a,b.c) 


By  making  use  of  definitions  (7.12a,b,c)  and  after  lengthy 
calculation  shown  In  Appendix  D  we  obtain: 

0^  V  )  ^  ^ 

'  xyyx' 


^XX  ll-V „.,v 


1 

fhB 

^1 

2^"lsin  f  1 

[^13 

|w,|  2 

A  } 

"'ll  1 


l*',l  iwjl  ,  ,v 

“  ~r^  Ti  iT5  ■  nr-  ~n  TT5  ^  \-U 


^  |wx|  B 


Iv^l  ® 


{7.13a,b.c) 


Extrapolating  the  results  found  from  equations  (7.11a,b,c, 
d)  the  stress  intensity  factors  k|j,  k^^y  can  be  com¬ 

puted  in  a  straight-forward  manner.  The  results  are  shown 
in  Figures  17-19.  Figure  17  shows  the  variation  of  kjj 
with  c/h2,  for  different  values  of  (h^/hg)  ratio,  kj^  in¬ 
creases  as  (h^/h2)  increases.  Figures  18  and  19  show  the 
variation  of  k^^^^  and  k^^y  with  respect  to  c/h2. 

7.3.2  Non-Sinqui ar  Behavior  at  the  Interface 

In  this  case,  the  characteristic  equation  (6.10)  has 
no  root  and  therefore  the  surface  displacement  derivatives 
are  bounded.  Since,  as  it  was  shown  in  Section  6,  the  dis¬ 
placement  derivatives  are  continuous  at  the  bimaterial  in¬ 
terfaces,  using  the  single-valuedness  condition  to  write 


-60- 


the  integrals  from  0  to  I,  equations  (6.2a,b)  take  the  form 
(see  [12]): 

If  71?  f  =  p(s)  0<s<t  (7.10) 

J  q  '  _1  ■'o 


where 


♦(t)  (0<t<h^,  0<r<h^l 
♦*(t)  (-h2<t<-c,  h,<» 


k(r.s) 


r  (1-v  V  ) 

I - P]  (*1 5  ( 0<x^  <h^ ,  0<s<h^ ) 

i  .  ♦  *  . 

I  •  *  *  (•h2^X2^“C» 

L  2Tri4^y 

r  * 

kl](x^,t)  (0<x^,  t<h^  ,  D<r,  s<h^) 

k*  fv  tl  if_L_+_!_l  °^s<h, 

12^1’  ^‘Mt-x,  t+x,J  hi<r<t 


k*  fv  tl  1  f-J—+  _lJ 
^l^^’^'"t-X2  t.X2  o<t<h,.0<r<h, 

>  ■’  I 

k22(*2’*^  (-h2<X2,  t<-c,  h^<r,  s<Jl) 


k,-j(Xi,t)  =>  k.j(x.,t)  -  k,j(x^.-t)  (i.j  =  1,2) 


Now  we  have  the  governing  equations  for  a  crack  imbedded 
in  a  non-homogeneous  material ,  obviously  with  power  singu¬ 
larity  -h  at  the  crack  tip.  Normalizing  equation  (7.10) 
by  means  of: 

r  »  tT  and  s  =  tic 

we  obtain: 

Go(T)dT  +  tj  kjj{T,ic)Gjj(T)dT  =  P|j(k)  0<ic<1  (7.11) 


i[’j- 


T+K 


Equation  (7.11)  reduces  to  a  set  of  linear  equations 
by  using  the  method  of  collocations  (see  [11]): 

1 


N 


1 


+  irtk„(T.,K:J 


Hg(Tj)  =  2Np^(<.)  i  =  1....N 

(7.12) 


where 


'*1 

Go(t) 


»  cos 


Holll 

/I  -T^ 


j  =  1,...N 
i  =  1,...N 


The  N  unknowns  H  (t-)  can  be  found  from  equation  (7.12)  in 

U 

a  straight-forward  manner. 


Defining  the  stress  intensity  factor  at  the  crack 
tip  as: 
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oi^ 


^(X2-i)  02ytX2’°^ 

we  obtain: 

2Y,-E„^r 

fcfa  -  -  .--  t  v-  H^O)  (7.13) 

Using  the  same  material  for  both  strips,  we  recover  again 
the  results  of  colinear  cracks  in  homogeneous  medium.  Fig¬ 
ure  20  shows  the  variation  of  kj^  with  c/h2.  k^^  increases 

as 


8,  COHCLUSIONS 

The  fracture  problem  of  layered  orthotropic  composite 
plates  has  been  studied.  The  following  results  have  been 
obtained: 

1)  Depending  on  the  elastic  constants,  orthotrOpic 

materials  can  be  classified  in  two  groups:  materials  of 

type  I  and  materials  of  type  II.  (A  different  formulation 
is  needed  for  each  combination.) 

2)  The  colinear  crack  solution  is  the  same  for  homo¬ 
geneous  isotropic  and  orthotropic  materials. 

3)  In  the  case  of  an  internal  crack  in  the  first 
layer,  the  stress  intensity  factor  k^  can  be  reduced  sig¬ 
nificantly  by  a  proper  selection  of  the  elastic  constants. 

4)  For  the  case  of  broken  laminates  there  is  a  sing¬ 
ularity  power  which  can  be  found  from  equation  (5.7).  The 
singularity  power  y  varies  between  0  and  1  for  different 
material  combinations. 

5)  For  a  crack  crossing  the  interface,  the  singular 
behavior  at  the  interface  disappears  for  some  material  com¬ 
binations.  In  this  case  the  crack  surface  displacement 
derivatives  are  bounded  and  continuous,  and  all  stresses 
are  bounded  at  the  bimaterial  interfaces. 
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9. 


RECOMMENDATIONS 


In  the  present  work,  a  general  formulation  of  the 
fracture  problem  of  layered  orthbtropic  composites  with 
periodic  cracks  is  given.  The  formulation  is  done  only 
for  the  case  where  both  materials  are  of  type  I.  Follow¬ 
ing  the  same  procedure,  the  problem  can  also  be  studied 
for  orthotropic  materials  of  type  II,  or  for  the  combina¬ 
tion  of  type  I  and  type  II.  The  dependence  of  the  singular 
behavior  at  the  interface  on  the  elastic  constants  can  also 
be  Investigated. 

In  our  formulation  the  thickness  of  the  adhesive 
bonding  the  layers  has  been  neglected,  the  §tudy  of  the 
effect  of  the  adhesive  also  can  be  recommended. 

A  more  realistic  approach  also  would  be  to  study  the 
problem  of  finite  number  of  strips.  But  this  problem  re¬ 
quires  lengthy  algebra. 

There  are  many  other  problems  to  be  studied  in  the 
fracture  of  composites.  We  hope  that  our  work  will  have 
a  small  contribution  in  the  study  of  these  problems. 
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FIGURES 
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Figure  8.  The  stress  intensity  factor  kj  for  the 
case  in  which  both  strips  contain  cracks 
(Combination  I), 


Figure  9.  The  stress  Intensity  factor  for  the  case 
in  which  both  strips  contain  cracks 
(Combination  I). 
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(Combination  II}. 
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laminate  Is  broken  (Combination  I  and  II). 


gure  i5.  The  stress  Intensity  factor  k,  when  the  firs 
laminate  Is  broken  and  the  second  contains  a 
crack  (Combination  I  and  II). 
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APPENDIX  A 

Definitions  of  the  material  constants: 

A  superscript  *  will  be  used  for  the  material  in  the  second 


strip.  The  constants  6.,  (i  = 

1,..,10)  and  Wj  (j  =  1,..,4) 

are  given  by  equations  (2.15), 

(3.5),  and  (3.8). 

,  '’yx''l®9 

Y,  =  1  +  ■  '  - 

Ti2  *  - 

•'Bs 

''vx‘'3®10 

T,  =■  1  +  ■  -- 

Yi3  =  sign(wi)B9  -  sign(w3)6io 

^3  *  “1  *  ''yx®7 

^14  7  >■  Y  Y 

^4  *  “3  *  V^S 

Y  =v 

5  vf7 

E  1  •*  V  V 

X  T  “x  xy  yx 

n  E  ,  *  * 

x  1  -  V  V 

V  '  '^xy  yx 

A,  =■  ^ 

'w 

Y  .V 

®  *y  V67 

^3  ^  V3  ■  ^7^4 

y?  "  ''xy’^l  ^  ®7 

*  ilr 

^4  “  ^iy3^7  ■  ®7^3 

^8  “  V3  ^  ^8 

^5  ^  ^iy4®7  ■  V3 

Y9  =  -1 

-  ^10 

^10  “  ^  ^s'^s 

* 

^7  ”  Yg  •  Yg^2 

V,i  =  - -^  +  sign(w,)Bg 
"  yf7  '  ^ 

* 

^8  °  ^9  ■  ^10^2 
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^9  “  ^3^7 

^24 

*  ^10^16  ■  V16 

^10  ”  ^3^8 

^25 

“  ^10^17 

“  ^3^9 

^26  ' 

*  ^12^16  '^10^18 

X,2  =•  ^6^4 

hi ' 

*  V20  ■  ^3^16 

^13  *  ^5^6 

hs '' 

=  XgX^g 

^14  “  ^6^3 

^29  ' 

*  ^12^19  '^13^17 

^15  ’  hh 

^30  ' 

'^13^18  '^12^20  ' 

^16  *  ^3^7 

^31  ' 

*  ^n^i6  ■  Vl6 

^17  “  ^3®7  ■*'  ^4®8 

^32 

,  X  X 
'  n  17 

Ai8  =  ^487 

^33 

“  Wl2  ■  ^11^18 

* 

^19  '  ^3®8  *  ^5^8 

^34 

=  XgXgg  -  Xi4X^g 

^20  "  ^5®7 

^35 

=•  XgXgl 

^21  “  ^3  ”  ^3®8 

^36 

“  ■  ^14^17  ■  ^21^ 

^22  “  ^3^7 

^37 

“  ^14^18  ■  ^22^12 

^23  ''  ^7^8 

^38 

2 

“  ^15^16  ■  ^7^9 
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^9 


X  X 
9  23 


X  a  X  X  .XX 
40  15  17  12  23 


^4  “  -  X,oX23 

^5  ’  ^13^23  '^15^19 


\l  “  ®7^12  "  ^15^18  ° 

*56“ 

■  V20 

^42  “  Vl6 

*57“ 

>3^19 

^3  “  ^3  ^17 

*58  “ 

-  X^X^g  (Xg  -  Yg) 

^44  “  Via 

*59“ 

■  VI6  ^Yg  -  X7) 

^45  *  ^10^16  "  ^11^16 

*60  “ 

Yg  (^4^ig  ■  ^5^17^ 

^46  “  ^llV  "  ^13^16 

*61  “ 

V7  <>10  -  ‘9) 

^47  “  "  ^11^19 

*62  “ 

^8  ^  Vl  7  ■'■  ^4^21  ^ 

* 

^48  “  ^14^16  ■  ^10^22 

*63“ 

-  *9  (Vs  ^ 

^9  *  ^13V  ■  ^14^20  °  ° 

*64  = 

-  ^7  (*10  -  *9) 

^50  *  ■  ^10^21 

*65  = 

*8  ^*4*23  *  ®7^7^ 

V  ”  ^14^19  *  ^13^21 

*66  = 

*7  {®7*19  -  *5*23) 

V  “  ^10^7  ■  ^15^16 

*67  = 

*  *4*16 

2 

*53  ^  ^15*20  ■  *13®7  ^  ° 

*68  = 

X  X 

5  16 
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r 


^69  “  Vl7  ■  Vl9 

^1  =  ■*■  ®8^5^^9  "^7^^ 

hz  “  ®7^^11  Wl3^  *  ^8^^3^12"^n^4^ 

^73  “  ^21^^10'^9^ 

^74  “  ^12^3  "  ^loH  * 

^75  “  ^9^5  ■  ^13^3  *  ^8^  V3‘^4^3^ 

^76  “  ^23^^10*^9^ 

^77  “  ®7^^10®7’^15^3^ 

^78  ^  ®8^^15V^^9^ 

A79  =  X2(X^63-X587) 

-80  “  ^24  *  hs  *  he  *  h?  ~  hs  *  h9 

^81  “  ■  ^70  ■  ^71  *  hz  *  S’9n(wi)B9(X73  +  Xj^  +  X^g) 

*  ~l^  ^he*  hi*  hz*  '  jw^  ^^43'^57'^^79^ 
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^82  “  ■  ^58  ■  ^59  ■  ^60  si3n(w^  )Bg(Xg:j  +  Xgg  +  Xgg) 

*  ^^64"^  ^es'^^ee^  ‘  |‘J^| 

^83  “  ^^70''’^71‘''^72^  '  7^  ®’9"^’^3^^10^^73^^74*^75^ 

12  12 

(^76‘^^77‘^^78^  *  '  |w  |  ’  ^^43‘^57'^^79^ 

12  3  3* 

^84  “  Y^  (^58'''^59''’^60^  ‘  ®’9'’^*'3^^10^^6l‘''^62‘^^63^ 

12  12 

Y^fw^  ^^64‘^^65'^^66^  *  "IwT'^^67'^W^69^ 

12  *  3  3 

X  -  -VsT 

85  Yi3Yi4Xgg 

X  .  ^8^82 

86  YT3Y^4Agg 

X  s  ^7^83 

■'^13‘'^14^80 

X  =  ^8^ 

^13^14^80 

^89  °  ■  ^70  "  Nl  ■  ^72  *  sign(w^)8gC^73  +  ^74  +  ^75) 


Y*  Y* 

■  1^  ^^76'^^77'^^78^  ■  ^2  |^*|*  ^^43"^57'^^79^ 

■  ^58  ■  ^59  "  ^60  ■  '*’  ^62  "'’  ^63^ 


*  '  ^  * 

■  ^^64'^^65‘^^66^  ■  h  ^^67'^^68‘^^69^ 
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Y  Y 

It  12 

^  ^^76'^^77'^^78^  ■*■  ‘iw*|  *  ^^43‘^57‘^^79^ 

Y*  Y* 

1*  12 

y*  y*  y*\<^ 

*  Y^  ■f;;|f  ^  V'*'W^66^  *’  *|w*l  *  f^67'*’W^69^ 

=  ^8^0 
■<^3^14^ 

J7*91 

■>'13^14^0 

/_V92_ 

’'13^14^ 

-  ^45  -  ^4g  *  ^47  +  slgnCw^Bg  (X^g  +  Xgg  +  Xg^) 

*  ^  1^  ^^52'^^54^^55^  ■  |w  I  ^^42'^^56'^^57^ 

*  I*  '  1 

■  ^1  ■  ^32  ■  ^33  '*’  "*■  ^35  ■*■  ^36^ 


■*■  ^^38'*’W^40^  '  }wj|'  ^  ’  ^^2  ■  ^43  *’  ^44^ 

1 

^  t^45  *  ^46  *  ^47>  •  ®’9"(”3^^10^^48  ^  ^50  ■"  ^51  ’ 

It  12 


Y^Tw^  ^^52'*’V'*’^55^  *  'itrr  ^^42  ■^^56‘^^57) 
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7^  (^31'*’^32'*’^33^  ■  ®’9"^“3^^10^^34‘*'^35''’^36^ 

12  12 


7^^  (^38'*'^39*^40^  *  Iwjl  ^42'^43‘^^44^ 

^97^7 

^13^14^80 

^1 3^14^80 

^99''^* 

^13^14^80 

.  ^100^8 
WT?80 

'  *  ^45  ■  ^46  ■  hi  '  S’9n(wi)89  (^43  *  ^50  *'  ^51^ 

Y*  Y* 

'  ^^52‘*’^54*^55^  ■  ^2  "j’*j''  ^^42*^56'*’^57^ 

’  ”  ^31  ■  ^32  ■  ^33'®’®"^*'1^®9  ^^34^^35'^^36^ 

■  ^^38'^^j9*^40^  ■  h  ’^^■^42’^43‘^^44^ 

’  ^  ^  #  s1gn(w3)6*o(X48n5o+X5^) 

II  12 


Y..*C 


Y*  Y* 

^1  7^  ~  '*' '  (^c9'‘'^i:A'*’^ct;)  ■*■ 


Y*  X 
II  2  P5 


55'  1^ 


(^42'^^56'^^57^ 


Y  *  *  * 

7^  ^^31*^32'^^33^  *  ^  sign(w3)S^Q(X3^+X35+X3g) 

12  12 

*  ^I’^T^'j^  (^38'^^39‘^^40^  *  ^|wfr^  ^■^42'^43^^44^ 
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:  ^105^* 
■'^13^14^80 

^106''^8 

^13^14^80 

*107^7 

^13^14^80 

* 

*108^8 

^13^14^80 


^113  “ 


.  UliIl^L 

^  «‘109 

/67  ^ 

“s 

^93^101 

Iliii  JwlL  X  x„ 

+  iwiiJwlL 

^  8® 

''85111 

.i!!iiJ[Mx  X 
>n5-7i^Wm- 

.i^iwtLx  X 

X  aliiLKLx  x„„  + 
8® 

+  i!liij«!lx  X 

.anAia^x  X 

✓B7  ^ 

-^^W,03 


.  l!iii  JilllL  X  X 
^ 

mKix  X 

IMJiLx  X 

iMi^X  X 

'  /8I  ^ 
JliiJ^lLx  X 


.iw^lwlix  X„1 

^  vif  ’ 

-?#■■■” 
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,  .JiUiiilj  1  ♦liiilil.x  X 

JT,  ^  ^  ® 

.liiLJiii,  ,  .jMJrlii  1 
^  ^  /5I 

x,„  ,  JlxlliiLx  X  -i^M-ljilx  X 

'»'  /S7  /5?  /S?  ^  ™ 
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APPENDIX  B 


Expressions  of  the  functions  used  In  equations  (4.14)  and 
(4.15): 


2Yi3«  ib 


-  s1gn(w3)6ioe'^*’2”'’^“’^^^”’^  'J)*(n)dn 

^12 


^11  n-(hi  -n)a*fs/  |w  I 


-  -^T-n-e 


3>  (ii(n)dn 


.  JS.  e-( 


2Yi3«ib[lwJ 

iLlA  ())*(n)dn 

y!,  Iwtl  J 


0  (a)  a  -  „  .1 .  f®C-(hl-n)a.^/|wi|.  Ill  g-(hTn)tx>/BT/|w3| 

'  2^3“  Jg  Yt2  i|>(n)dn 

-  —1—  e'(^2-n)a>^/|wTl.  Ip.  e-('’2-n)a>^/’wti  ^*(r,)dn 

Zy^jO  ib  ^12 

R2{a)  =  2^  f  s1gn(wi)B,e-(>'l-'’5“’^/l«*l 

-  p-  s1gn(w3)6:oe^^'’l-^’“'^/l‘'3l  $(n)dn 
^12 

-  _4_  f‘’rsign(w*)B:e-(>'2-n)«-^/lwtl 


.  ^  4,(^)dn 

f'’|Ve-(»'2-n)a^/|w:i 

2Y*3a  lb  Iwtl 

-  -L-e  •^^2-'i)«-^/|wtj  ♦*{n)dn 


f(o)  =  X24f^(fli)  +  X25f2(a)  +  ^26^3^“^  *  ^27^4^“^ 

it  ■  It 

f|(a)  =  tanh{w^ah2)tanh(w20ih2) 
f2(a)  a  tanh(w2ah2)tanh(w^ah^) 
fjCa)  3  tanh{w2ah2)tanh(w20ih|) 
f^{a)  a  tanh(w*ah2)tanh{w2ahi) 

fg(a)  =  tanh(w*ah2)tanh(w^ah^ ) 

fg(o)  =  tanh(w^ah^ )tanh{w2ah|) 

g{a)  =  Xjitanh(w*ah2)  +  X22tanh(w^ah.j )  +  X22tanh{w2ah^ ) 

h{a)  =  +  XjsfgCa)  + 

ni(a)  “  X2Q‘f^(a)  +  ^3g‘^g'.^)  ^40^6^*^^ 


n(a)  =  -  X^2taiih(w^ah2)  -  i^jtanhCw^ah^ )  +  X^^tanh(w2ahi ) 
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go(a)  =  A^gtanh(w2ah2)  +  X^gtanh{w2oh^)  +  X^7tanh(w^oh^ ) 

"0^“^  =  X^2^anh(w2ah2)  +  AgglanhCw^oh-l )  +  Xgytanh{w^ah^ ) 
g.|{a)  =  Xggtanh(w2ah2)  +  XggtanhCw^ahg)  +  Xggtanh(w^ah^ ) 

h^(a)  =  Xg^f^(a)  +  Xg2f2(a)  +  Xggfg(a) 


ni,{a)  =  Xg4f,(a)  +  Xggf2(a)  +  Xggf5(a) 
it  it 

n^(a)  =  Xg7tanh(w2ah2)  +  Xggtanh(w^cih2)  +  Xggtanh(W|ah^ ) 

92(0)  a  Xygtanh{w3ah2)  +  X^^tanh(w*ah2)  +  Xy2tanh{w2ah^ ) 

^  ^74^3^®^  *  ^75^4*®^ 

in2(a)  =  Xyg'F^  (ot)  +  Xy^f2(ot)  +  Xygf^(ct) 

n2{a)  =  X^gtanh(wgah2)  -  Xgytanh(w*ah2)  +  X^gtanh(w3ahi ) 
Expressions  of  the  functions  k.{j=l,8)  used  in  Eqs.  (4.17): 

V 

,  jcosh(w,ax, )  f 

‘'l^^l’®^  *  2Y,3f(a)[cosh(w^ah,)  ^7  -92(“)  +  sign(w,  )Bgh2(a) 


Y..*C 


COShCWjQX^ ) 


+  ■  "tw"!"  "2^“^  coshlwgah^}  ^8  '^l 

y 

+  sign(wi )  Bgh,  (a)  +  "l 
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Y,  ,  ' 


+  sign{w.j)ggh(a)  +  i/^ m(a) - — n(oi) 


kg{x2.o)  - 


y  cosh(ii/^ax2)  * 


'^Yi3TI«J  kosh(w^oh2) 


^177^7  9o(“) 


Y  Y 
11  ^  *11 

Y  "  Y 
12  12 


sign(w3)B^ghjj{a) 


Y  Y  Y  ^  COSh(w~ax,)  *  Y 

■ '^^TTT  ^  "~7  >  Yg  7^g(a) 

*12  '"s'  '3'  cosh(w2ah2)  ^12 


Y..  Y. 


Y..«C 


-  ^  sign{w3)$^gh(a)  -  ni(a)  +  "(“) 

1  fcosh{wrax-)  *  f  *  * 

’^^"2*“^  ”  7T77tLc7  *  .7  ^7  -9of“)  -  sign{w,)e9h^(a) 


2yi  3^(01)  j^osh(w^ah2) 

Y*  Y* 

/I  °  2  iw!i 


’  Tvr  -  ^2  “PT  "0^“^  -'7'  *  ^  ,  ^8  -9(“) 

IwJ  IwJ  cosh(w3ah2) 


cosh(w3ax2)  * 


1 

C0Sh(w*aX2) 

2y*  f(a) 
13 

cosh(w*ah2) 

-  sign(w*)Bgh(o)  -  ni(a)  -  n(a) 

iwjt  lw,l 

,  cosh(w^ax-)  *f  Y* 


-  ^  sign(w*3)g;,  h„(a)  .  ^  A  "0^“^ 

'12  Yjj  l"3l 

Y*  ^  ^  COSh(w,ax,)  *  Y*  Y*  * 

+  5-  n  (a)  + - — - Yg  g(a)  +  signCwj) 

|w  1  ^  cosh(w*ah2)  ®  y*z  Y*  ^ 

3  J  0  4^  — I** 

*  Y*  Y*  Y*  ^ 

B,jjh(a)  +  X,  ^  ^  in(a)  +  — L  n(a) 

’  -  ’  ‘  Y*  Iw!!  iM 


The  kernels  K.(j  =  l,4)  used  in  Eqs.  (6.15a,b): 

J 
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IC-|{n.y)  =  57—  r  -  y-  ^  e'lwslay/'^  cosah^sincirida 

‘  ^^13  ^o[  '  ‘ 

+  j^|*  -  y3g2(a)  +  Y49i  (a)  +  sign(w^  )eg  Y3h2{a) 

]  Y,  f  1  Y  vr  f 

+  Y4hi  (a)  +  1 — r  Y3in2(<=‘)  *  Y4m^  (a)  -  ■•■■—p  Y3n2(a)  +  Y4n^  (a) 

IWjl  |Wj| 

f  f  ' 

4,  g-(hi-n)a>^/|w3l  ^  Y3g2{a)+Y4gi(a)  - ^  sign(w3)3,Q  Y3h2{a) 


Y4h3(q)  Y3ni2(a)+Y4niT{a)  +  -|^  ' 


Y3n2(a)+ 


K2{n,y) 


-(h2-ri)a^/|w*|_  g-(h2+r))ci>^/|wi|  I 


Y3g2(a) +  Y4gi (a)  -  sign(Wy)Bg  Y3h2(a) +  Y4hi (a) 


♦  r 

T  ^  ■ 

-  >^-V-  Y,nip(a) 

IM 

*  Y4'^ (a)  -  *- 4-—  Y3n2{a) +Y4n^(a)  +'  ^'^3 

]  j  i  . 

.  e-(h2+n)a/^/|wtl  ,  .  L 

Y*  .  J  J 

+  sign(w3)B^g  Y3h2(a) 
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Y*  Y*^  T  Y*  f 

+  Y4h,(a)  +  ^^-\^Y3m2(a)+Y4'ni{a)  +  YjngW 

12  i’"3l  |W,| 


f 

ICjCn.y)  =  — |°°  Y-ji  sinah-iSinanda 

*  C  -[Ygtanh(w^ah^)g2(a) 


Y^Qtanh(w2ah^)g^(a) 


s1gn(w^  )6g  Ygtanh{w^ah^ )h2(a)  +  y^o*'®'’*’^'*3“*’i 


+  -j-i|-  |^Ygtanh(w^ah^)m2(a) +Y^Qtanh(w2oh^)ni^(a) 

Y  ^ 

-  Ygtanh(w^ah^)n2(a) 

+  Yiotanli(w3ah^)n^(a)  Ygtanh(w^ahi)g2{a) 

J)  iH 


+  Y|o<^3nh(w3ah^)g^(a)  -  sign(w3)8^Q  Ygtanh(w^Bh^ )h2(a) 


YiQtanh{w3ahi) 


Y, ,  Y. 


h^(a)  -  /g7  Ygtanh(w^ah^)m2(a) 


}VK 


|w,| 


Ygtanh(w^ah^  jPgCoi)  +  Y^Qtanh(w2cih^  )n^  (a)l 


sina 


TW. 


iiS^da 


rr  -*  "i 

-  f 

1-  Ygtanh(wiah^)g2(a)  +  Yio<^3nh(w3ah^)g^(o)|  -  sigh(w*)6g 


Y9tanh(wiah^)h2(a) +Y^Qtanh(w3ah^)h^(a)j  -  ^YgtanMw^cih^ ) 

]  Y*  X  C 

"^(a) +yiota''f’(w3“h^)'"l(“)j  -  ~[^[~  ^YgtanhCw-ah^ )n2(a) 

+  Y^gtanh(w3ah^)n^(a)||  + 


a'(h2+ri)a^/| 


Y^gtanh(w3oh^)g^(a)| 


Ygtanh(w^ah3)g2(a) 


y 

*  ^  sign(w3)6*g|Ygtanh(w^ah^)h2(a) +Y^Qtanh(w3ah3)h^(a) 


V*  *A 

j.  VT  I  ^2^1 

’'e « w 

*  ^12  |wj 


Ygtanh(w^ah^  )nt2(a)  +  Y^gtanh(w3ah^  )m^  (a) 


Y*  v'b*  ^ 

+  lU-Es—?. 


Iwtl 


Ygtanh{w^ah^  )n2(a)  + 


SI  nay 
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APPENDIX  D 


Derivation  of  the  Stress  Intensity  Factors 


D.  1  Case  of  Internal  Cracks: 


The  stress  intensity  factors  are  defined  as: 


=»  Tim  v^(x^-a)  (x^.O) 


X;|-*a 


k|j  “  lim  /2(x2-b)  02  (xg.O) 


X2-*'b 


From  Eq.  (4.16a)  we  can  write; 

y  V  ) 
xy  yx' 

where  a°y(x^,0)  is  a  bounded  function. 


1^  +  a°y(x,  .0) 


,iii/2 


*(t).-ELtL  =  _Fitk 

(t-a)*5(t+a)** 


(D.1a,b) 


(0.2) 


Define  the  sectionally  holomorphic  function, 

r 

From  [14]  we  obtain: 


<.(z!  .  1  p  ^  « 


t|»(z) 


F(-a)e 


Tri/2 


(2a)’*(z+a)’*  (2a)’*(z-a)*^  ° 

Using  (0.3)  to  evaluate  (D.2)  and  with  the  definitions 
given  in  (D.la.b)  we  have: 


(D.3) 


-no- 


XF(a)  =  - 

''xy^yx^  ''xy^yx^ 


and  similarly 


t  i f*(i,) .  -  F^C) 

b  (T_v  V  1  /b  fl-v  V  )  “ 

xy  yx^  '  xy  yx' 

D;2  Case  of  Broken  Laminates 


From  Eq.  (5.2b) 


2YTaE!  fh 


where  O2y(x2«0)  *  bounded  function. 


♦(t)  = 


(h|-t*)‘'^  (t-h^)^(t+h^)^ 


Define: 


♦('I  ■  H '  W  « 

-h] 


F(-h,)e"^  1  F{h,)  T 

^(2)= - J - ! - ! - ! — _+i{){z)  (D.e) 

(2h-j)^sinirir  (z+h^)^  (2h^)^sinwY  (z-h^)^ 


Using  (D.4),  (0.5)  and  (0.6)  we  obtain: 


(h^)\*^Fo(1)  r  N.lZ/gT 
^’■''xyV^®’"^  -  jwtlKI  ^ 


Y  *  ^102 


|wtl|wi| 


-m- 


*  -  X  t 

hslIwTi 

Y  n04 

Iwsllwtl 

Y 

- 

k|j  has  the  same  expre  rion  as  the  one  we  derived  in  the 
previous  case. 


D. 3  Case  of  a  Crack  Crossing  the  Interface 


D.3-a  Singular  Behavior  at  tiie  Interface; 

In  this  case  the  definitions  of  the  stress  intensity 
factors  are: 


•‘xx  * 


and 


“xy"  ;4JV  y^xy('’l'y> 

From  Eq.  (6.2b)  we  can  write: 

rh,  ,*, 


"  (l.v*  V*  ):r  ^  ^  <^2yU2-0) 

'  xy  yx' 

where  a2y(x2.0)  is  a  bounded  function. 


(D. 7a ,b.c} 


(D.8) 


(h2-c)®sintr/2  (z-c)** 


_ 1  ^  i/)*(z) 

(h2-c)*^simTB  (z-h2)®  ° 


Using  (0.9)  and  (D.8),  Eq.  (D.7a)  yields  to: 

*  * 


Ic  -  F*(c)>^ 


(0.9) 


2^-^;4E*(h2-c)'s  . 

I^k  * - T-V^— 

^  (l-v  V  )  ° 

'  xy  yx' 

Separating  the  singular  parts  of  Eqs.  {6.15a,b)  we  have: 


/■hi 

aix(hi,y)=  j  K^5{n.y)<|)(n)dn 


+ 1  ■*■ 
^Ixy^^’r^^  "  K3s(Ti.y)<fr(n)dn 

+  I  %3(n,y)<!i*(n)<ln  +  T°j^y(y) 
where  o°j^(y)  and  T°j^y(y)  are  bounded. 


(D.lOa.b) 


then 


(2h^)®sinirB  (2+h^)®  (2h^ )®simTB(z-h, ° 


(D.n) 


Using  (D.ll),  (0.9),  (D.IOa.b)  and  the  definitions  (0.7b,c) 
after  lengthy  algebra  we  obtain: 


and 
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The  derivation  is  similar  to  the  one  done  for  the  case  of 
internal  cracks.  Thus,  we  have; 
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